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1 Introduction 



The semileptonic and radiative decays of heavy mesons play an outstanding 
role for determination of the parameters of standard model, in particular the 
quark mixing parameters. 

The upcomming and currently operating B factories BaBar at SLAC, 
Belle at KEK, LHCB at CERN and B-TEV at Fermilab as well as the 
planned r-Charm factory CLEO at Cornell make precision tests of Stan- 
dard Model (SM) and beyound SM ever more promising [1]. Especially, a 
strigent test on the unitarity of the Cabibo-Kobayashi-Maskawa (CKM) mix- 
ing matrix [2] in the SM will be made by these facilities. Accurate analysis 
of exclusive semileptonic B decays are thus strongly demanded for such pre- 
cision tests. Beauty decays proceeding through the flavour changing neutral 
current (FCNC) transition have been identified as a valuable source of infor- 
mation on standard model and its extensions [3] (for a review and complete 
set of refrences see [4, 5, 6]). One of the physics programs at the B factories 
is the exclusive B decays induced by FCNC transitions. These decays are 
forbidden at the tree level in the SM and occurs at the lowest order only 
through one-loop (penguin) diagrams. The experimental tests of exclusive 
decays are much easier than those of inclusive ones. 

The theoretical understanding of weak decays of hadrons, and the mea- 
surements of the corresponding CKM matrix elements are consistently ham- 
pered by the presence of long distance QCD effects that are responsible for 
the binding of quarks into hadrons. These effects are hard to evaluate in 
a model independent way, and so tend to bring large uncertainities to the 
theoretical predictions for the weak decay amplitudes. They appear in the 
calculation of the matrix elements of the weak Hamiltonian operators, be- 
tween the initial and final hadronic states. 

The dynamical content of hadronic current matrix elements is described 
by Lorentz-invariant form factors. The calculation of these form factors re- 
quires a non-perturbative treatment and are source of large theoretical un- 
certainities. Knowledge of these form factors is essential for the description 
of semileptonic and non-leptonic weak decay processes and in particular for 
the experimental determination of CKM matrix elements [2]. Many theo- 
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retical tools have been applied for the determination of form factors, i.e. 
quark models, QCD sum rule, and lattice calculations. The machinery of 
sum rules, lattice calculations has the advantage to be directly based on the 
QCD Lagrangian, but is quite involved. Quark models, on the other hand 
are less directly connected with the QCD Lagrangian, but give a vivid pic- 
ture of what is going on and allow an easy application to different processes 
and at quite different kinematical regions. However, models presented so far 
lacked full relativistic covariance with respect to quark spins, and it is the 
spin structure which determines the ratio of different form factors we are 
concerned. The knowledge of these ratios is of paramount importance for 
the analysis of experimentally measured decays. So, the relations between 
different form factors that will hold under certain conditions or approxima- 
tions can then be very useful. One of the motivations of the present work 
is to make use of Ward Identities to relate various form factors in a model 
independent way. Such relations will reduce the number of uncertain quanti- 
ties, and improve the accuracy of the theoretical predictions. Moreover, they 
may help us understand better the general features of the underlying long 
distance QCD effects. 

As mentioned earlier, the theoretical understanding of exclusive decays 
is complicated mainly due to the non-perturbative hadronic form factors 
entered in the long distance nonperturbative contributions. These calcu- 
lations of hadronic form factors for semileptonic B decays have been in- 
vestigated by various theoretical approches. Form factors based on lattice 
calculations [7, 8], and on light cone sum rules [9, 10, 11], currently have 
uncertainties in the 15% to 20% range. A variety of quark model calcula- 
tions [12, 13, 14, 15, 16, 17] exist. A number of other approaches [18, 19], 
such as dispersive bounds and experimentally constrained models based on 
Heavy Quark Symmetry (HQS) also exists. In the standard model, the short 
distance contributions are dominated by the top-quark and long distance 
contributions by form factors. 

The aim of the present work, as mentioned earlier, is to relate the vari- 
ous form factors in a model independent way through Ward identities. This 
enables us to make a clean separation between non-pole and pole type contri- 
butions, the q 2 — > behaviour of the former is known in terms of a universal 
function £j_ (q 2 ) = g + (q 2 ) introduced in the large energy effective theory 
(LEET) of heavy (B) to light (p) form factors [20]. The corrections of order 
(m/M) 2 , where m is the light meson mass and M is the heavy meson mass, 
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in LEET were obtained in [21]. The LEET approach was first proposed by 
Dugan and Grinstein [22] and further developed by Aglietti [23]. 

The residue of the pole is then determined in a self consistent way in terms 
of g + (0) or £j_ (0) [a sort of duality]. Thus our approach has a predictive 
power as the residue of the pole will give information about the couplings of 
B* (1 _ ) and B* A (1 + ) with Bp channel. The form factors are then determined 
in terms of g + (0), A (0) (for which much information is already known) and 
pole masses M B * : M B * A (where M B * is known and estimates exists for M B * A ), 
and M'/M [where W is the radial excitation of M (B* and B\)\. This ratio 
is 1.14 (1.12) [24]. 

The information about the hadronic form factors thus obtained are used 
to calculate the decay rates and the CKM matrix element V u b- Because \ V u b\, 
the smallest element in the CKM mixing matrix, plays a crucial role in the 
examination of the unitarity constraints and the fundamental questions on 
which the constraints can bear [25]. To extract \V u b\ from exclusive channel, 
the form factors for the channel must be known. The form factor normal- 
ization dominates the uncertainty on \V u b\. The q 2 dependence of the form 
factors, which is needed to determine the experimental efficiency, also con- 
tribute to the uncertainty, but at much reduced level. There have been two 
exclusive \V u b\ analyses by the CLEO Collaboration: a simulatneous mea- 
surement of the B — > nW and the B — > plv transitions [26], and second 
measurement of the B — > plv rate [27]. The branching fractions obtained 
were 

B (b° -> p~l + v) = (2.5 ± 0At° 5 7 ± 0.5) x 1(T 4 [26] 

B (B° -> p~l + v) = (2.69 ± 0.411°;*; ± 0.5) x 1(T 4 [27]. 

The results of the two analyses are largely statistically independent, and they 
have been combined 

\V ub \ = (3.25 ± 0.14+°;^ ± 0.55) x 10" 3 , 

where the errors arise from statistical, experimental systematic, and form 
factor uncertainities, respectively. A recent study by BaBar [28], obtained 
the following results for the branching ratios and \V u b\ and 

B (b° -> p~l + v) = (3.29 ± 0.42 ± 0.47 ± 0.60) x 10^ 4 , 
IK&I = (3.64 ± 0.22 ± 0.25+°;^) x 10" 3 , 
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where the uncertainities are statistical, systematic, and theoretical respec- 
tively; and CLEO [29], 



where the errors in CLEO analyses are statistical, experimental systematic, 
theoretical systematicbased on lattice QCD and LCSR uncertainities, and 
plv form factor shape uncertainly, respectively. 

The future for the exclusive determinations of |14&| appears promising. 
For both lattice and B factories, B — > nlu appears to be a golden mode 
for future precise determination of |V^| while the B — > plv mode will be 
more problematic for high precision because the broad width introduces both 
experimental and theoretical difficulties. 

This paper is organized as follows: In section 2, we give the B (p) — > 
V (k, e) Iv current matrix elements. In section 3, we discuss the Ward Identi- 
ties and develop the relations between form factors which result in to reduce 
the number of unknown quantities. In section 4, pole contribution of vari- 
ous form factors are discussed and relations among different decay constants 
with the help of Ward Identities are obtained. In section 5, form factors are 
obtained which are necesary for the calculation of decay width and branching 
ratio. In section 6, we discuss the decay distribution in terms of helicity am- 
plitudes and the extraction of \y u b\ from the experimental data. In section 7, 
numerical analyses is made and our predictions regarding various observables 
are compared with other approaches. Section 8, summarizes our conclusions 
and predictions. 



We are interested in the exclusive semileptonic decays of pseudoscalar meson 
(B) into vector meson (p). The structure of hadronic current in semilep- 
tonic decay must be constructed from the available four-vectors, which are 
momenta and spin-polarization vectors. The Lorentz-vector or axial vec- 
tor quantities thus formed have Lorentz invariant coefficients (form factors) 
that are functions of q 2 = (p — k) 2 . The polarization vector of the meson 
p (k, e) leads to a hadronic current with three form factors (in the limit of 
zero charged lepton mass). 




2 Current Matrix Elements 
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In the case of B (p) — > p (k, e) Iv decay, there are seven form factors for 
pseudoscalar to vector transitions, i.e. V, ^-0,1,2, and i*i,2,3 [4, 5, 12] 

2ie " Ua " VWfiVtf) (1) 



(V(k,e)\V»\B(p)) 
(V(k,e)\A»\B(p)) 



Mb + My 
(M B + M v )e*»A 1 (q 2 
e* ■ q 



Mb + My 



-2M, 



£ ■ q 



(k+pfA 2 (g 2 ) 
A 3 (q 2 ) ~ A (q 2 ) 



where V** = wy^b, A^ = wy^b, and 

/ 2 N Mb + M v 
A3 (q ) = 

with A 3 (0) = A Q (0), 



2Mi 



v 



M B - M v 



2M* 



v 



A 2 (g 2 ) 



(2) 



(3) 



(4) 



(V (k, e) \uia^q v b\ B (p)> = -ie^elKp^ (q 2 ) 
V(k,e)\ma^ 1 5 q u b\B(p)) = [(m 2 b - M 2 ) e*» - (e* ■ q) (k + pf] F 2 (g 2 ) 



<f - 



M 2 -M 2 



withFi(O) =2F 2 (0). 

Again the term proportional to q^ in Eqs. (2-5) only plays an important 
role for the case / = r and is not relevent for the rate in case of produced 
lepton pair I = p, e. The same is true for the form factor F 3 . However, it 
does play a considerable role in the Ward identity. 



3 Ward Identities 



The Ward Identities, which are derived in appendix A, relate various form 
factors [30] 

(V (k,e)\uia^q u b\B(p)} = - (m b + m q ) (V (k,e) \u^b\ B (p)) (6) 
[V (k, e) \uia^ 5 qi/ b\ B(p)) = (m b - m q ) (v (k, e) \u^b\ B (p)) 

+ (p + kf (V (k, e) u-fb B (p)\ (7) 



-iq^{V(k,e) 



B(p) 



i (m b + m q ) (V (k, e) vrfb B (p) ) (8) 
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Using the Ward Identity (6) in Eqs. (1) and (4), and comparing the coeffi- 
cient, we obtain 



m b + m q 



■2V(q 



_. , (9) 

M B + My ^ > K } 

Again, using the Ward Identity (7) in Eqs. (2) and (5), and comparing 
the coefficients of e M * and that of [we do not compare the coefficients of 
(p + kY as heavy quark approximation has been used in obtaining the second 
term of Eq. (6b)] we obtain 



F 2 (q 2 ) = 



m b - m q 
M B - M v 



F,(q 2 ) = -(m b -m q ) 2 -^f{A 3 (q 2 )-A (q 2 ) 



(10) 
(11) 



The results given in Eqs. (10), (11) and (9) are model independent be- 
cause these are derived by using Ward Identities. 

In order to obtain a universal normalization of the above form factors at 
q 2 = 0, we define 

(V (k, e) \iua a pb\ B (p)) = -ie afipa e* p [(p + kf g + + q a gJ [ - q ■ e* (k) e afipa (p + kf q a h 

-i [q a £/3 P ar£* p ip + k f q T - a <-> (3} hi 

-i [(p + k) a e PpaT e*P (p + k) a q T -a^ 0] h 2 (12) 

Since 



"a/3, 



(13) 



we can write 

(V(k,e) \iua^ lb b\B{p)) 

Using 



and 



__c.fival3 



e^P(V(p>)\iua a pb\B(p)) (14) 







e 0123 




9s 


9 a p 


9a 


Aivafi _ _ 






& 
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< 


9a 



-1 
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h 



we obtain from Eq. (12) 

(V(k,e)\tua^q ul5 b\B(p)) = e^[(M 2 B - M 2 ) g + + q 2 g 

-q-e*(k) [(p + kf g + + q»g_] 
+q-e*(k) [q 2 (p + kf - (M 2 B - M 2 ) <f 

-e** 4 (M 2 B -M 2 )\ 2 

+£*" (2M 2 B + 2M 2 - q 2 ) q 2 h 2 

+q ■ e* (k) (p + kf (M 2 -M 2 - q 2 ) h 2 

+q ■ e* (k) q» (-M% - 3M 2 + q 2 ) h 2 (15) 

Equating Eqs. (5) and (15) and comparing various coefficients 

F 3 (V) = -g_-(M 2 -M 2 )h-(M 2 +3M 2 -q 2 )h 2 



(1 



F 2 (q 2 ) = g+ -q 2 h-(M 2 B -M 2 -q 2 )h 2 



M 2 -M 2 



F 3 (q 2 ) 



= 9+ 



M 2 B -M 2 



9- ~ 



Q 



Ml-Ml- 

B v M 2 B -Ml 



(2M 2 B + 2M 2 - q 2 ) 



(1 



Further, from Eq. (12) 



V (k, e) \iua a ^b\ B (p)) = -ie aPpa e*^q^2 [g + - q 2 h x - (m 2 b - M 2 ) h 2 

(18) 

Comparing Eqs. (18) and (4), we obtain 



F 1 (q 2 ) =2[g + - q 2 h x - (M 2 B - M 2 ) h 2 



(19) 



Note that above relations ensure that F 2 (0) = (1/2) Fi (0). Equation (10), 
by using the value of F 2 (q 2 ), becomes 



M« 2 ) = 



M B -M v 
m b - m q 



9+ + 



M 2 -Ml 9 - 



Ml -Ml : r 

B V Mg — My 



(2M 2 B + 2M 2 - q 2 ) 



h 2 \ (20) 
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Equation (19), by using the value of F\ (q 2 ) given in Eq. (9), becomes 

M B + M v 



V (f) = 



m b + m q 
Rearranging Eq. (3), we get 



9+ ~ q 2 h - (Ml - M 2 V ) h 2 ] 



(21) 



A 2 (q 2 ) 



M B - M v 
M B + M v 
M B - My 
M B + M v 



mM*) — 



2Mi 



v 



M> 

2My 



V 



M B - My 



A 3 (q 2 ) 
As (q 2 ) - A (q 2 ) 



2M, 



v 



M B - My 



Ao (q 2 ) 



g + - q 2 h 
m b -m q y 

2M V A f 2 x 

-m^ 4 (« ) 



M 2 B 



M 2 + 



Ml -My 



(3Mf + 5M V - 2q 2 ) 



(22) 



where we have used Eqs. (11) and (16). 

The merit of above relations is that the normalization of the form factors 
Ai, V, Fi and F 2 at q 2 = is determined by a single constant [g + (0) — {M B — M v ) h 2 (0)] 
while those of A 2 and A are determined by [g + (0) — (M| — My) h 2 (0)] and 
■In ('0). 



4 Pole contribution 

Only hi, <?_, h, and A get pole contributions from B* (1~), B\(l + ), and 
B (0~) mesons. These are given by 



hi 



pole 



9— \pole 



h\ 



Ao (q 2 ) 



pole 



pole 



1 9B*Bp- 

2 M 2 B , 

9B* A Bp- 



B* 
T 



R 



V 



1 



1 - q 2 /M 2 B * M%. 1 - q 2 /M 2 B « 
•t 



£ A 



R S A 



M 2 B * 



1 - q */M 2 



B* 
T 



Ml, 



q 2 IM\ 
l 



M 2 



9bb p - 
M n 



h 



1 - g 2 /M 
q 2 /M 2 



M 2 



q 2 /M 2 



— Rq- 



q 2 /M 2 



l-q 2 /M 2 "l-q 2 /M 2 



(23) 
(24) 

(25) 
(26) 
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On the other hand, only g + , g_, and Aq get contribution from quark 
A-graph. Therefore we shall put h 2 = in what follows. 
The above coupling constants are defined as follows: 



B* + (rj,q)p-(e,k) £° (p) ) = ig B *B P -e apiia e* a qW» 

(0 \iua^ v b\ B* (77, q)) = {q^ v - q v r]^) (27) 
B* A + ( V ,q)p-(e,k)\B°(p)) = ig^-e* ■ 77* - f B ^ Bp -q ■ e*k ■ 77* 

(0 \iua^ v b\ B* (77, q)) = f T A e l _ lual3 r] a q fJ 

Now if in the Ward Identity (6), we take the matrix elements between (0| 
and \B*), we obtain 

(0 \%ua^q v b\ B* (77, q)) = - (m b + m q ) f B .rf (28) 

where (0 \iu^^b\ B* (77, q)) = f B *rf. Thus using Eq. (27), we obtain 



fB* m b + m q m b + m M B 



It is easy to see that same relation is obtained if one uses heavy quark spin 
symmetry. Similarly, if in the Ward Identity (7), we take the matrix elements 
between (0| and \B* A ), we obtain on using that (0 |m75&| B* A (77, q)) = 0, 

(0 \iua^q vlb b\ B\ (77, q)) = (m b - m q ) f B *y, (30) 

where (0 |iu7 /i 75&| B* A (77, q)) = f B * A rf. By using the identity (13), we obtain 

B * A _ m b -m q 

fr --^rr^\ ( 31 ) 

A 

5 Form factors and determination of coupling 
constants 

The form factors V (q 2 ) [Eq. (21)], A 1 (q 2 ) [Eq. (20)], and A 2 (q 2 ) [Eq. (22)] 
can be written as 

v/ m b + m q L v/ v / J 
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M B + M V \ f 2 x n g 2 



M B " v M 2 *l-q 2 /M B * 

M B -M v j / 2 x g 2 (~ ( 2\, ( 2\ 



W-^^T3OT ^2) 



A, (g 2 ) = ^^{. + (.V^K^)}-M^^(g 2 ) 

M B + My \ , 2 x D g 2 1 1 2My / - 

" M^My ( g ) ~ Ra Mf A TZ qV Ml« A ) - M^My^ W 

The behaviour of g + (g 2 ), g_ (g 2 ) and A (g 2 ) near g 2 — > is known from 
LEET 

The pole terms in the relations (32), (33) and (34) are expected to dom- 
inate near g 2 = M|» or M|. . Therefore, we cannot expect these relations, 
obtained from Ward identities, to hold for all g 2 for which we use the param- 
etarization, suggested by above behavior, near g 2 = and near the pole 

F ( n 2\ = F ( Q ) (on) 

V y ) (l-g 2 /M 2 )(l-g 2 /M' 2 )' 1 ' 

where M 2 is M B * or M|* and M' is the radial excitation of M. This param- 
eterization takes into account potential corrections to single pole dominance, 
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presumably arising from radial excitations of M [30], as suggested by dis- 
persion relations [31, 32] (a detailed discussion on this point one can find in 
[30]). It also takes care of off-mass-shell-ness of couplings of B* or B* A with 
Bp channel [30] by replacing the pole contribution 



R 



with 



R(l- M 2 /M' 2 ) 



l-q 2 /M 2 """"" (1 - q 2 /M 2 ) (1 - q 2 /M' 2 ) 
Since g + (q 2 ) and g_ (g 2 ) do not have a pole at q 2 = M B , , we get 

M B + My 



(40) 



This becomes 



q 2 =M 2 B * 



M f 



1 

T 



Rv = —-gB*B P - fr 



-^9B*Bp-fB 



-Rv ■ 



(41) 



Similarly, using parameterization (39) for the first term on the right hand 
side of Eq. (34), 



R 



D 



1 , .b* a _ 1 , M B -M V 

--^jB* A Bp-JT - -^3B A Bp- ^72 JB\ - 



9+ (0) 



1-M 2 /M' 2 



(42) 



We cannot use the parameterization (39) for A 1 (q 2 ) since near q 2 = 0, A\ (q 2 ) 
behaves as g + (q 2 ) [1 — q 2 / (M B — M v )]. This suggests the following 



9 + (0) 



1 - 



(l - q 2 /M 2 B1 ) (l - q 2 /M' 2 A ) \ M 2 - M 2 J 
so that, since g + and g_ do not have pole at q 2 = M|» , 

1-M 2 /(M 2 -M 2 ) 



(43) 



A 1 (q 2 ) [l - g 2 /M^ 



q 2 =Ml, 

B A 



9+ (0) 



1 - M£. /Mg, 
1 



2 __ M B -M V f 
M 2 -M 2 B *~ 9B * Bp - Mh * A ~M\ — M 2 



v 



Most of the estimates of g + (0) suggest that [33] 

g+(0) =0.29 ±0.04 [11] 



(44) 



(45) 
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Table 1: B- mesons masses in GeV [12, 24] 







M 


M'/M 


M B 


cr 


5.28 


1.14 


M B * 


\- 


5.33 


1.14 


M b * 


i+ 


5.71 


1.12 



For numerical work we shall use the S-meson masses given in Table 1. 
Thus we have the prediction from Eq. (41) 

2 (1.26 ±0.17) 
9b*b p - = t (46) 

JB 

This prediction [/# = 180 MeV] is consistant with its various estimates e.g. 
QCD sum rules give it \[2 (11) GeV" 1 with gs*B P = V^9b*b p {> (where gs*B P ° 
has been estimated in [32]) and the constituent quark model (CQM) predicts 
it to be [31] 

^{K^ + i)f^}^ (9 - 75)G ^" (47) 

The predictions (42) and (44) cannot be tested at present since we do not 
have any information on fs* A - However, the ratio of S-wave and D-wave 
couplings is predicted to be 

9b* a b p - M^-iMl-M*) 



fs* A B P - 2 

g B * A B P - = 2.66 x f B * A B P - GeV 2 (48) 
The different values of F (0)'s are 

V(0) = ^^, + (0) (49) 

A 1 (0) = g + (0) (50) 
A , x M B + M V , . 2M V . . 

M0) = wt^ 9+{0) -M^k M0) (51) 
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where g + (0) is given by (45). As far as A (0) is concerned we can get some 
infor mation by writing 



A (q 2 ) = A (q 2 ) + A (q 2 ) 



pole 



a ( 2\ u <f/M 2 B {l- M 2 /M' 2 ) 

~ A ° ^ ) +R \i- q yM B )(i- q yM>£) (52) 

Using the parameterization for A (q 2 ) similar to (39) and the same argument 
as before [cf. Eq. (41)], we get 

/ M 2 \ 

A)(0) = -lo (0) /.'„ 1 ' 11 ] 



9_BBp_ f ( _M] 
M n 13 { M' 2 



Using the p-dominance 



so that 



M 2 1 

9BBp- = j^j^ (54) 

^^'■frNw (55) 



where F p (0) takes care off mass shell-ness of p-meson and is given by [31] 



/ 9pirnfp \ 



-1 



Fp{0) = [viw =a80 ' (56) 



exp 



where we have used 

I 9pnirfp 



1.22 ±0.03. (57) 



Using / B = 180 MeV, f p -/M p = 198 MeV, 

A (0) = 0.26. (58) 

Although LEET does not give any relationship between £y (0) and £j_ (0) but 
due to some numerical coincidence in the LCSR expressions for £y and £j_ 
[20, 38] 

£|| (0) * a (0) = (0) (59) 
14 



so that, from Eq. (37) 



A (0) = 1.100+ (0) ~ 0.320 ± 0.044 (60) 

consistent with the above estimate (58). Another estimate of A Q (0) in light- 
cone sum rules (LCSR) approach is 0.372. 

We summarise the form factors which we shall use for our numerical work 
[M v = M p = 0.770 GeV] 

v ( q2 ) = (i - eiMD (i - <?/m§) (61) 

A ^ = (l - g 2 /M%) (l - 7JM^7) i 1 ~ Ml - M v ) (62) 

A ( Q 2) = MO) 2M^ A (0) 

2{q) (l - ^/JWl.) (l - 52/M«.) M B -M v {l- q ^/Ml){l- q ^M'l 

where 

V (0) = 0.332 ± 0.046, 

Ai(0) = 0.29 ±0.04, (64) 

A 2 (0) = 0.389 ± 0.054. 

These form factors (5) are plotted in Fig. 1, where the comparison of the form 
factor vs. lattice data [7, 43, 44] and calculations within different approaches 
is given. Solid line-our result, dotted line-LCSR [11], dashed line-quark 
model [14], dash-dotted line-lattice constrained parametriation of [7]. 



6 Decay distribution 

Four independent kinematical variables completely describe the semileptonic 
decay B — > pW, where the vector meson decays to two pseudoscalar mesons, 
p — > 7T7T. The four variables commonly used are q 2 = (p — k) 2 and the three 
angles (0i, 9 P , and x)- The angle 0i is measured in the W* (or IV) rest frame, 
it is the polar angle between the charged lepton and the direction opposite 
to that of vector meson. The angle 9 p is measured in the rest frame of vector 
meson, where the pseudoscalar mesons are back to back, and it is the polar 
angle between one of these mesons. The angle x is chosen to be azimuthal 
angle between the plane of flight of lepton and of pseudoscalar mesons. 
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The differential decay rate for B — > plv, p — > txtx can be expressed in 
terms of these four variables g 2 , 9i, 9 P , and x [45, 46] 



dT (B — > plv,p — > 7T7r) 
<ig 2 d cos # p d cos 9i dx 



8(4tt 



\ V u b\ 



Ml 



-B(p 



7171) 



cos 9if sin 2 (g 2 ) 



+ (1 + cos ^) 2 sin 2 ^ 



+ 4 sin 2 9i cos 2 61 



—4 sin 0/ (1 — cos ^) sin 9 p cos # p cos xH+ (g 2 ) H (q 2 ^J 
+4 sin 6^ (1 + cos 9i) sin9 p cos9 p cosxH- (g 2 ) H (g 2 ) 
-2 sin 2 9 X sin 2 9 p cos 2x#+ (g 2 ) H (g 2 ) } 



where 



1 



(M 2 B + Mj-q^y 
AM 2 B 



M 2 



(66) 



is the magnitude of three-momentum of p in the rest frame of B and is a 
function of g 2 . 

Bescause the parent meson has spin zero, the vector meson and the W* 
must have the same helicity. The amplitudes for the helicities 0, +1, and -1 
are proportional to H (g 2 ), H + (g 2 ), and if_ (g 2 ). The detailed dynamics of 
the hadronic current are described by the variation of these helicity ampli- 
tudes with g 2 . Equation (65) incorporates the V — A structure of the leptonic 
current, as well as the assumption that the mass of the charged lepton can 
be neglected. 

The helicity amplitudes can in turn be related to the two axial form 
factors, Ai (g 2 ) and A 2 (g 2 ), and the vector form factor V (g 2 ), which appear 
in the hadronic current (2) and (1) respectively, 



H (g 2 ) 



2M pV ^ 



(Mi - g 2 - M 2 ) (M B + M p ) A l (g 2 ) 



4M 2 



M B + M, 



and 



H±(q 2 ) = (M B + M p )A 1 (q 2 ) T 



2Mb \K\ 
M B + M. 



-A 2 (g 2 ) 
(67) 

(68) 
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The form factors A\ (q 2 ), A2 {q 2 ), and V (q 2 ) are dimensionless and relatively 
real, since CP is conserved in these decays and there are no final state strong 
interactions [16]. We note that, while A 2 contributes only to H and V 
contributes only to H±, Ai contributes to all three helicity amplitudes. At 
high q 2 (small |k p |), each of the helicity amplitudes is dominated by A\. At 
q 2 = 0, the over all factor of q 2 in (65) causes all the contributions to the 
decay rate to vanish except that from \H (q 2 )\ 2 , since H contains a factor of 
l/Vq 2 and at q 2 = q 2 .^ = (M B — M p ) 2 , the over all factor of |k p | vanishes, 
causing the rate to do likewise. 

Integrating (65) over the angles 9 P and Xi we obtain 



d 2 T 
dq 2 dE l 



G 2 F I v u b I q 2 

128vr 3 M 2 B 



{(l + cos^) 2 |#-(V) 



+ (1 



cos^) 2 \H+ (g 2 )f + 2 (l - cos 2 ^) \H (q 2 )\ 2 \ (69) 



where 61 is the polar angle between the p and the lepton / in the (lv{) center 
of mass system, and 



2M B \K\ cos °i = M b~ M p +<f ~ 4M B E l 



(70) 



The Ei (or cos 6{) integration of Eq. (69) can be done trivially and results 
in the differential g 2 -distribution 



dT_ 

dq^ 



G 2 F \V„ 



■ub 



96tt 3 



K\ ^f^P - ttt) (^)l 2 + |^ + + \Ho {q 2 )\ 2 } 



(71) 

By numerically integrating Eq. (69) over lepton energy Ei in the interval 
< Ei < 2.6 and over q 2 in the interval 



< q 2 < 2Ei 



'M 2 B 



M 2 
p 



2M B Ei 



M B - 2Ei 



the width of B — > plu turns out to be 



r (B -> plu) = 15.82 ± 2.18 \V ub \ 2 ps-\ 



(72) 



where we have used the form factors given in Eqs. (5). 
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The branching ratio of B° — > p l + v is measured by CLEO collaboration 
[27], B expt (B° -> p~l + v) = (2.57 ± 0.29t° il ± 0.4l) x 10" 4 , where the errors 
are statistical, systematic, and theoretical. The CKM matrix element \V u b\ 
can be obtained from the branching ratio B expt (B° — > p~l + v) using 



V ub \ = 



= 3.24 x 10~ 3 . (73) 



We use r B o = 1.548 ± 0.032 ps [47]. Recently, the BaBar Collaboration [28] 
predicted the branching ratio B expt (B° -> p~l + v) = (3.29 ± 0.42 ± 0.47 ± 0.60) x 
10~ 4 and |Kb| = (3-64 ± 0.22 ± 0.25±g;|§) x 10~ 3 , where the uncertainities 
are statistical, systematic, and theoretical respectively. By using the branch- 
ing ratio predicted by BaBar, and decay width (72), we obtain 

\V ub \ = 3.67 x 10~ 3 . (74) 

Another recent study by CLEO [29], gave the branching ratio B (B° — > p~l + v) = 
(2.17±0.34i|];^±0.41±0.01)xl0- 4 , |K&| = (3.00 d= 0.21±g;|§ ±g;^ d= 0.28) x 
10~ 3 , where the errors in CLEO analyses are statistical, experimental system- 
atic, theoretical systematicbased on lattice QCD and LCSR uncertainities, 
and plv form factor shape uncertainity, respectively. By using this branching 
ratio given by CLEO [29], we obtain 

|Kb| = 2.98 x 10~ 3 . 

In this study of B — > plv decay, there is only one free parameter i.e. g + (0), 
and we have used the value of this predicted by a LCSR calculation [11] given 
in Eq. (45). The error appear in our calculation from this parameter. The 
error shown in our result of T (B — > plv) in (72) is only 13.8%. 



7 Numerical Analyses and Discussions 

The values of the helicity amplitudes and form factors at q 2 = and q 2 = g max 
are given in Tables 2, 3, and 4. The variation of form factors at q 2 = and 
q 2 = g max is shown in Table 4. None of the helicity amplitudes changes sign 
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Table 2: A comparison between different calculations is given between the 
helicity amplitudes H + (GeV) and their ratio H-/H+ at q 2 = 0. 





H + (0) 


#-(0) 


H_ (0) /H + (0) 


V¥H (0) 


Our 


0.256 


3.253 


12.714 


8.737 


QM [14] 


0.163 


2.942 


18.014 


8.202 


LCSR [11] 


0.055 


3.103 


56.767 


10.158 


Lattice [7] 


0.0898 


4.224 


47.020 


10.914 



Table 3: A comparison between different calculations is given between the 
helicity amplitudes if_, H+ (GeV) at q 2 = g^. 





H+ (<Zmax) 


H (?max) 


Ho (<?max) 


V¥H (g^ ax ) 


Our 


2.362 


2.362 


2.362 


10.652 


QM [14] 


2.952 


2.952 


2.952 


13.313 


LCSR [11] 


2.782 


2.782 


2.782 


12.549 


Lattice [7] 


2.826 


2.826 


2.826 


12.744 



over the entire physical range of q 2 as shown in Fig. 2. Although H + and 
Ho both get contributions from two form factors which enters with opposite 
signs. The sqaure of the helicity amplitudes is plotted in Fig. 3. 

First, H- is generally dominent over H + , while the condition H = H + 
is forced by kinematics at q 2 = f^axi these two amplitudes go their own way 
and differ by more than a factor of 12.7 at q 2 = 0, a comparison over this 
point is given in Table 2. Again kinematics forces H = H_ = H + at 
q 2 = g^ax' but the amplitudes soon separate so that in the region near q 2 = 
the amplitude Hq, which contain a factor of l/y/q 1 , very much dominates the 
others. This feature is independent of modest changes in the form factors. 

We have used vector-current form factor pole at 5.33 GeV and axial- 
vector current form factor pole at 5.71 GeV. Note that the poles are very 
close to the edge of physical region. This means that on the one hand, they 
should indeed dominate the behavior of the form factors near q 2 = g^ax- The 
change in the q 2 range in different form factors give quite distinct results as 
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Table 4: Variation of form factors at q 2 = (q 2 = g^) is shown. A com- 
parison of the ratio of respective form factors at q 2 = g^ ax over at q 2 = is 
presented. 





V 


A 1 


A 2 


no) 


Ai(0) 


A 2 (0) 


Our 


0.332 (2.605) 


0.29 (0.390) 


0.279 (1.136) 


7.840 


1.346 


4.06 


QM [14] 


0.308 (1.915) 


0.257 (0.488) 


0.242 (0.956) 


6.218 


1.901 


3.955 


LCSR [11] 


0.338 (2.011) 


0.261 (0.460) 


0.223 (0.818) 


5.948 


3.135 


2.062 


Lattice [7] 


0.458 (2.102) 


0.356 (0.467) 


0.342 (0.898) 


4.587 


1.310 


2.628 



shown in Table 4. 

The H + (H) helicity amplitudes gets destructive (constructive) contri- 
butions from the A\ and V form factors. The near zero behavior in H + 
arises, since V with its coefficient 2Mb |k p | / (Mb + M p ) in Eq. (68) causes 
its contribution to this halicity amplitude to quickly catch up to that of A\ as 
we move away from q 2 = g^ax- Of course this increases if_ correspondingly 
according to Eq. (68). 

A similar situation pertains to the helicity amplitude H , where Ai and 
A 2 destructively interfere, as shown in Eq. (67). This pulls the whole longi- 
tudnal portion of the decay rate (Tl = T ) down and it is far smaller than 
the transverse one (T T = T + + r_), see Table 5. This agrees with the obser- 
vations of Gilman and Singleton [45], and Korner and Schuler [16] (as shown 
in Eq. (19) of Ref. [16]) but contradicts the recent observation of Ali and 
Safir [48] in the LEET approach. 

As an exercise, we check the relation among the form factors which holds 
true in large energy limit (LEL) (for a recent discussion see Refs. [20, 49, 50, 
51]), namely 

Ai (q 2 ) = 2E V M B = M 2 +M 2 - q 2 
V(q 2 ) (Mb + My) 2 (M B + M v f ' 

In Fig. 4, we plot the ratio of B — > p form factors A\ (q 2 ) /V (q 2 ) and 
compare them with the LEL [20] Eq. (75) (dashed line), with the result of 
Ref. [39, 21] (dotted line) and the lattice results [52], whereas the solid line 
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Table 5: Different decay rate values are calculated and compared with dif- 
ferent approaches. The comparison of Yl/Yt is presented and compare with 
other approaches. 





r_ 


r+ 


r L = r 




^Total 


Tl/Tt 


Our 


9.662 


0.059 


6.095 


9.721 


15.82 


0.63 


QM [14] 










15.8±2.3 


0.88±0.08 


LCSR [9] 










13.5±4.0 


0.52±0.08 


Lattice [7] 










16.533 


0.80^ 


ISGW2 [13] 










14.2 


0.3 



Table 6: The values of the V/A\ at different q 2 is given and compared with 
HQET-LEET approach and with recent Lattice data. Also ratio T + /r_ is 
given. 



Observable 


V/A, 
q 2 = 


V/A, 
q 2 = 14.53 GeV 2 


V/A 1 
q 2 = 14.68 GeV 2 


r+/r_ 


Our 


1.146 


2.96 


3.00 


0.006 


HQET-LEET [20] 


1.29 


2.63 


2.66 


0.02 


Lattice [52] 




2.05 


2.02 
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is the result of this calculation. The knowledge of this ratio has important 
consequences: on the one hand, its measurment at q 2 = will provide a test 
of various approaches; on the other hand it provides the ratio r + /T_ of the 
width to helicity eigenstates A = ±1. In Table 6, the values of the V/A 1 
at different q 2 is given and compared with HQET-LEET [20] approach and 
with recent Lattice data [52]. Also ratio r + /T_ is given in Table 6. 



8 Conclusions 

We have studied B — > plui decay by using Ward Identities. The form fac- 
tors have been calculated and found that their normalization is essentially 
determined by a single constant g + (0). Then we use a parametrazation [Eq. 
(39)] which take into account potential corrections to single pole dominance 
arising from radial excitations of M (where M = M B * or M B *J, as suggested 
by dispersion relations [31, 32]. This also takes care of off mass shell-ness of 
couplings B* or B* A with Bp channel [30] as indicated in Eq. (40). We took 
the value of g + (0) = 0.29 ± 0.04 from a recent LCSR calculations [11]. And 
predict the value of g B * Bp = 2 (1.26 ± 0.17) / f B GeV" 1 [cf. Eq. (46)], which 
agrees very well with CQM calculation (47) and the value obtained by QCD 
sum rules a/2 (11) GeV^ 1 . Also, we have predicted a relation between S- 
wave and D-wave couplings g B ^ Bp - = 2.66 x f B ^ Bp - GeV 2 given in Eq. (48). 
The estimate for A (0) ~ 0.320 ± 0.044 is consistant with the value given 
in Eq. (58) and with the LCSR value 0.372 [11]. We summarize the form 
factors in Eq. (5) and their values at q 2 = in Eq. (64). 

A detailed analysis of form factors and helicity amplitudes is presented. 
Total decay rate is obtained which is 15.82 \V u b\ 2 ps _1 [cf. Eq. 72], which 
results in predicting the CKM matrix element \V u b\ 

\V ub \ = 3.24 x 10" 3 

by using CLEO branching ratio result [27], and 

\V ub \ = 3.67 x 10~ 3 

by using BaBar branching ratio result [28]. While using the branching ratio 
of recent study of CLEO [29], we obtain the value of \V u b\ 

\V ub \ = 2.98 x 10~ 3 . 
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Our results of \V ub \ agree very well with the respective results of \V u b\ of 
CLEO [27], and the recent results of BaBar [28], and CLEO [29]. The error 
shown in our result of T (B — > plu) in (72) is only 13.8%. 



A Derivation of Ward Identities 

We start with the relation 

(V (k, e) \iua^q vl5 b\ B (p)) e'^'^ = - (V (k, e) \d v (u (x) a^ l5 b (x))\ B (p)> 

(76) 

We can replace d v by the covariant derivative D v to take into account the 
strong interaction of the the quark field. We can write the right hand side of 
Eq. (76) as 

- (V (k, e) \D v u (x) [-i-f-f + 75& (x)\B (p))-{V (k, e) \u (x) [i 7 " 7 " ~ lM (x)\ B (p)) 

(77) 

Using Dirac equation 



/Db (x) = —imiyb (x) , u (x) p — im u u (x) 

the relation (77) becomes 

(m b - m u ) (V (k, e) \u (x) 7^756 (x) \B[p)) 
-i(V(k,e) \D»u (x) 75ft (x) \ B (p)) 
+i (V (k, e) \u (x) l5 D»b (x) \ B (p)) 
= (m b - m u ) (V (k, e) \u^ lb b\ B (p)> e"** 
-i{V{k,e)\D» (u (x) l5 b(x))\B(p)) 
+2i (V (k, e) \u (x) -f 5 D% (x) \B(p)) 

Using now the linear momentum commutation relation 

>*, O (x)] = -iD^O (x) , P£ = j d 3 x : g f (x) D^q (x) 



and 



ly(k,e)\P» 
P" [5 (p)> 



fc"(y(A;,e)| 
p"|fl(p)), 



(78) 



(79) 



(80) 
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the last two terms in Eq. (79) become 

-<f (V (k, e) \u l5 b\ B (p)> e"** + 2 (V e) |u 75 6p£| 5 (p)) e"^ (81) 

where in the last term we have used that P£ \ V (k, e)) = as V e) does 
not contain the quark 6. In the heavy quark effective theory mj is taken 
to infinity in a way that p b /m is held fixed, but the four momentum of the 
light degree of freedom are neglected compared with nib [53]. This enables 
us to identify p^ with p M . Hence from Eqs. (76), (79) and (81), we obtain 
the Ward identity [Eq. (7)] [2p - q = p -> k] 

(V(k,e)\iua^q ul5 b\B(p)) = (m b - m u ) (V (k,e)\u>f>y 5 b\B (p)) 

+ (p^ + ^) (y(A;,e) |u 7 5&|£(p)) 

(82) 

Similarly, 

(y (fc, e) [iua^q^bl B (p)) = - (m b - m u ) (V (k,e) \u^b\ B (p)) 

+ (p^ + F) (F(A;,e)|fZ6| J B(p)) 

(83) 

However, the matrix element (V (k, e) B (p)) vanishes due to parity con- 
sideration. Thus, we obtain the Ward Identity (6) given in the text. 
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B Figure Captions 

1. Comparison of the form factor vs. lattice data [43, 44] and calcula- 
tions within different approaches is given. Solid line-our result, dotted 
line-LCSR [11], dashed line-quark model [14], dash-dotted line-lattice 
constrained parametriation of [7]. 

2. Comparison of Helicity amplitudes vs calculations within different ap- 
proaches is given. Line description is as in Fig. 1. 

3. dT/dq 2 distribution for each of the three terms in Eq. (71): (a) the 
terms proportional to \H_\ 2 and |i?+| 2 , (b) the terms proportional to 
\Hq\ . Line description is as in Fig. 1. 

4. Ratio of B — > p form factors, A 1 (q 2 ) /V (q 2 ), are plotted with lattice 
data points [52]. The solid line is our result. The dashed line is the 
LEET prediction [20] and the dotted line is the prediction of [39, 21]. 
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